Copyright © 1980 American Telephone and Telegraph Company 

The Bell System Technical Journal 

Vol. 59, No. 6, July-August 1980 

Printed in U.S.A. 



Propagation Characteristics of 
Double-Mode Fibers 

By L. G. COHEN, W. L. MAMMEL, C. LIN, and W. G. FRENCH 

(Manuscript received January 21 , 1980) 

We study propagation characteristics of double-mode fibers by 
measuring group delay differences of time-resolved discrete pulses of 
the lp (01) and LP (11) modes as a function of wavelength in the 1.06- 
to 1.6-\im region. The technique has been used to make the first 
measurements of zero intermodal dispersion in the zero material 
dispersion region of graded- index double-mode fiber guides that can 
have larger diameters than standard single-mode fibers. The results 
can also be used to deduce information about the fiber's refractive- 
index profile, such as A, the normalized peak refractive index differ- 
ence. 

I. INTRODUCTION 

This paper describes transmission characteristics of double-mode 
fiberguides which can have larger diameters than single-mode fibers 
without sacrificing bandwidth. Wide bandwidths are achieved in dou- 
ble-mode fibers because there is a wavelength at which the group 
delays are the same for their two propagating modes. 

The normalized frequency parameter, V, can be used to determine 
the number o f modes in a fiberguide. When V < V co = 
2.41 Vl + (2/g), then only one mode, 1 designated lp(01), 2 propagates 
in a guide with a graded refractive-in dex pro file described by n = 
no Vl — 2A(r/a)*. Since V = (2ira/\ ) V2re Are, it is evident that the 
core diameter, 2a, and the core-to-cladding index difference, Are = re(0) 
— re(a), may be varied over a wide range while still maintaining single- 
mode operation. 

However, when Vis significantly smaller than the cut-off value, V co , 
the field propagating in the fiber becomes loosely guided and trans- 
mission losses increase due to bending and cladding imperfections. To 
minimize material dispersion effects, bending losses, and splicing dif- 
ficulties, single-mode fibers should be designed with V a V,„ at a 
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wavelength as close as possible to a system carrier wavelength in the 
A = 1.3-pm region. 3 ' 4 This is achieved by designing core diameters to 
be about 15 /im and by having small index differences (An « 0.0015) 
between the core and cladding. Single-mode fibers designed for the 
0.8-jum region have somewhat smaller core diameters that are typically 
10 jLtm. 5 

Double-mode fibers could be used as an alternative to single-mode 
fibers if the group velocities of its two propagating modes were equal- 
ized at a system wavelength in the 1.3-jum region 6 or in the 0.8-/nm 
region. 7 Then intermodal dispersion would be zero in waveguides with 
larger diameters than single-mode fibers. However, practical systems 
applications would require precise double-mode profile control since 
the zero double-modal dispersion wavelength depends on the shape of 
the profile. 

The remainder of this paper describes a technique for determining 
the wavelength at which intermodal dispersion becomes zero in double- 
mode fibers. These are the first results reported in the 1.3-pm wave- 
length region where material dispersion effects are also very small. 
The technique can also be used to deduce the wavelength at which a 
single-mode fiber becomes double-moded as well as physical properties 
of the fiber such as its effective g and An. 

Measured results are compared to theoretical predictions based on 
exact numerical solutions for LP (01) and lp(11) propagating modes of 
the scalar differential wave equations for the electromagnetic fields. 
The numerical analysis includes the effect of dips in the refractive- 
index profile which occur at the center of the fiber core. 

II. THEORETICAL RESULTS 

Figure 1 shows intermodal propagation characteristics of double- 
mode fibers with two different power law refractive-index profiles. 
Normalized modal group delays t, defined by r = t (1 + At), are 
plotted versus V-number for values ranging from about 2 to 6. These 
results were obtained from exact numerical solutions for propagating 
modes of the scalar differential wave equations for the electromagnetic 
fields. The horizontal dashed lines in the figures correspond to the 
time delay, t /L, for a plane wave propagating in a medium with the 
peak fiber core refractive-index, n. The two solid curves in each figure 
represent delay characteristics for lp(01) and lp(11) modes. The single- 
mode limit for a fiber occurs at the cut-off V-number for the lp(11) 
mode which is V co = 2.41 for step-index fibers (Fig. la). Only two 
modes propagate for ^-numbers to the left of the dashed vertical line 
where the lp(21) mode becomes cut off. 

Zero intermodal dispersion occurs at the intersection point of the 
modal delay curves. Theoretical results, to be given in more detail in 
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Fig. 1 — Normalized group delays, t defined by t = to(1 + A ?), vs normalized 
frequency, V, for lp(01) and lp(11) modes in fibers having step (g = oo) and graded (g 
= 2.5) refractive-index profiles. The horizontal dashed lines in the figures correspond to 
the time delay, t /L, for a plane wave propagating in a medium with the peak core 
refractive-index, n. The vertical dashed lines indicate the cut-off V-number for the third 
propagating mode. 

separate papers, 8,9 show that this point is almost entirely dependent 
upon the fiber's refractive-index profile shape. The equalization 
V-number is 3 for step-index fibers. Figure lb applies to fibers with 
g = 2.5 profiles and shows the advantage of using graded-index double- 
mode fibers. The cut-off V-number is V co = 3.34 for single-mode 
operation and the double-mode regime corresponds to V< 5. However, 
in practice, double-mode operation can be extended to V-values 10 to 
20 percent higher than the third mode cut-off where the lp(21) mode 
is still very lossy due to the finite cladding thickness. Zero lp(01)- 
lp(11) intermodal dispersion in g = 2.5 fibers occurs when V = 5.5, 
which is almost two times larger than V co for single-mode operation 
and the equalization V-number for step-index fibers. This theoretical 
result implies that double-mode fibers with zero modal dispersion can 
be fabricated with core diameters that are twice as large as for single- 
mode fibers with step-index profiles. 

Fibers fabricated for this study were made by the modified chemical 
vapor deposition (mcvd) technique. 1 Their cores consisted of GeOi- 
doped silica and the cladding of pure silica. For a fiber with intentional 
index-grading, the core was deposited in several layers with different 
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Fig. 2 — Profiles of refractive-index vs radius, which are representative of fibers used 
in this study. The rms width, w, of the central index depression was about 25 percent of 
the core radius. Its depth, d, varied from about 50 to 100 percent of the maximum core 
radius. 



Ge0 2 concentrations. For the step-index fiber, the core composition 
was kept constant for each layer, but diffusion and/or dopant vapori- 
zation yielded some unintentional profile grading. Vaporization from 
the innermost layer caused refractive-index dips at the center of the 
core of each type of fiber. Figure 2 shows analytic profiles, with g = 
2.5, of refractive index versus radius which are representative of fibers 
used in this study. The normalized profiles in Fig. 2 are given by 
1 - N(r), 

r" + d-exp(-(r/w) 2 )(l - r M ) - b 



N(r) = 



1-6 



(1) 



where b is chosen so that the minimum of N(r) is zero. The rms width, 
w, of the central index depression was about 25 percent of the core 
radius. Its depth, d, varied from about 50 to 100 percent of the 
maximum core index. 
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Figure 3 shows results of numerical studies to determine the effects 
of central refractive-index dips within power law g profiles on the 
propagation characteristics of double-mode fibers. The group of three 
curves on the left show g plotted versus the lp(11) cut-off V-number 
for single-mode operation. The solid curve applies to profiles without 
dips. Note that g decreases from 8.0 to 2.5 as V co increases from 2.7 to 
3.3. The dashed curves show the effect of a 100-percent deep index dip 

whose width is 10 percent of the core radius and the ( ) curves 

show the effect of a 100-percent dip whose width is 25 percent of the 
core radius. These dips tend to increase the width of the propagating 
fields. A fiber without a dip having zero modal dispersion at V = 4 has 
a larger g than a fiber with a dip which has zero modal dispersion at 
V = 4. This makes the effective profile g seem larger than the 
fabricated g. 

The group of three curves on the right of Fig. 3 show how refractive- 
index dips affect double-mode operation by reducing the V-numbers 



g s 




Fig. 3 — Power law profile-^ vs V-number. The solid curves apply to profiles without 

dips. Data shown as ( ) apply to profiles with w = 1, d = 1. Data shown as ( ) 

apply to profiles with w = 0.25, d = 1. The group of three curves on the right are plotted 
vs the V-number for zero double-modal dispersion. 
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at which LP (01) and lp(11) modes are equalized. Note that the effect 
on double-mode operation is substantially greater than on single-mode 
operation. But again the dips make the effective profile g appear larger 
than the fabricated g values. For example, a 25-percent wide, 100- 
percent deep dip makes a fabricated profile g of 2.5 act like an 
undistorted profile with g = 5. This indicates that the effect of the 
profile dip could be partially compensated for by a fabricated profile 
with g < 2.5. 

When intermodal dispersion between LP (01) and lp(11) modes is 
zero, the fiber dispersion is due to intramodal material and waveguide 
effects for each of these modes. These effects are illustrated in Fig. 4, 
which is a plot of the dimensionless dispersion coefficient D versus 
V-number, for the parameters listed, where D is defined by the equa- 
tion 



dX c\ 



(2) 



The effects of material and waveguide dispersion are interrelated, and 
the total chromatic dispersion is the sum of the standard material 
dispersion common to all modes and waveguide dispersion that in- 
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Fig. 4 — Dimensionless dispersion coefficient D vs V for fibers with g = oo and g = 2.5. 
D can be converted to units of ps/Km-nm through (no/c\)b(8\)D. Material dispersion 

effects are illustrated by ( ) data for A = 0.001. The solid curves show waveguide 

dispersion effects for lp(01) and lp(11) modes. The fiber parameters a. A, Vo were 
calculated to make the zero double-modal dispersion wavelength, Ao, equal to the zero 
material dispersion wavelength, 1.282 jtm. 
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eludes material effects. In Fig. 4, the ( ) curve shows the material 

dispersion coefficient, D m versus V. The solid curves show the wave- 
guide dispersion coefficient, including material effects, for the lp(01) 
and lp(11) modes. The parameters are chosen so the zero material 
dispersion wavelength coincides with the zero intermodal dispersion 
wavelength. For germania-doped silica with A = 0.001, the material 
dispersion is zero at Ao = 1.282 /im. This leaves the core radius a to be 
chosen so the V-number for zero material dispersion coincides with 
the V-number for zero intermodal dispersion. For example, fibers with 
g = 2.5 profiles have zero intermodal and material dispersion at Ao = 
1.282 jum if their index difference is A = 0.001 within a core radius 
a = 17.6 nm. The dimensionless dispersion coefficient D can be 
converted to pulse-broadening in ps/Km-nm by 3.76D ps/Km-nm, 
through (lio/cA) A(5A)D, with no = 1.45, c = 3 X 10 5 Km/s, A = 1.282 
[im, and with (5A) . being the spectral bandwidth of the source in 
nautical miles. Therefore, waveguide dispersion effects, illustrated in 
Fig. 4, should cause less than 0.5 ps/Km-nm pulse broadening at the 
zero intermodal dispersion wavelength for double-mode fibers. 



III. EXPERIMENTAL RESULTS 

Time-domain pulse-delay measurements were made with the Uni- 
versal Fiber-Optic (ufo) measurement system described previously. 10 
The optical source provides subnanosecond pulses at 1.06 /im and at 
any wavelength within the 1.1- to 1.6-/mi spectral range. Transmission 
data are automatically acquired and processed by a microcomputer 
graphics terminal, which is now an integral part of the measurement 
system. 

Transmission measurement results in Fig. 5 show that LP (01) and 
lp(11) modes can be time-resolved and verify that zero modal disper- 
sion occurred in a double-mode fiber which was 500 meters long and 
had a graded refractive-index profile with an effective g ~ 3. Impulse 
responses are displayed on a 200-picosecond/div horizontal time scale 
for six wavelengths. The two modes in the fiber propagate as distinct 
pulses whose relative peak amplitudes could be changed by varying 
the input beam launching conditions. At A = 1.125 /im, the first (least 
delayed) peak represented power propagating in the LP (01) mode 
because it was associated with a nearly Gaussian far-field radiation 
pattern. The second peak was associated with the lp(11) mode which 
has a doubly peaked, far-field radiation pattern. LP (01) to lp(11) modal 
dispersion characteristics are determined by measuring the intermodal 
group delay differences between the pulse peaks. Note that, as the 
wavelength increases, the time delay between peaks decreases until 
the pulses coalesce into a single narrow pulse atg= 1.24 fim wavelength 
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Fig. 5 — Impulse responses are displayed on a 0.2 ns/div horizontal time scale for six 
wavelengths in a 0.5-Km-long graded-index fiber with an effective g ~ 3. Zero double- 
modal dispersion occurs at A = 1.24 /mi. The fiber is single-mode at A = 1.45 /im because 
the lp(11) mode is cut off. 



where the fiber V-number is V = (2iia/\) V2n An = 4.25. For wave- 
lengths longer than X = 1.24 fim, the lp(11) mode is delayed less than 
the lp(01) mode, and the modal pulses become further separated as 
the wavelength increases to 1.28 /xm. Finally, at A = 1.45 urn, the LP(ll) 
mode is cut off and the fiber is single-mode. Note that the double- 
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mode fiber response at A = 1.24 /zm is no more dispersive than the 
single-mode response at A = 1.45 jtim. 

Figure 6 summarizes results from several different double-mode 
fibers. The intermodal group delay difference between the lp(01) and 
lp(11) modes is plotted versus V-number. Zero intermodal dispersion 
occurs when the curves intersect the horizontal axis. The solid theo- 
retical delay curve, S, applies to step-index profiles with Arc = 0.016. 
The maximum lag time of the lp(11) mode behind the lp(01) mode is 
5t a 500 A/i ns/Km and can be used to estimate An in step-index 
double-mode fibers. This enabled us to deduce that A/i = 0.016 and 
0.010 in fabricated step-index fibers S-l and S-2, respectively. 
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Fig. 6— Intermodal time delay differences St = (t(01) - t(11)) vs V-number for 
double-mode fibers. Curve S shows calculated results for a step-index fiber. Curves S-l, 
S-2 are results for fabricated fibers with nearly step-index profiles. Curves G-l, G-2 

correspond to fibers with intentional index grading. The ( ) delay curves were 

calculated for g = 2.5 profiles which are nearly optimal for double-mode operation. 
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In general, the relative delay between lp(01) and lp(11) modes is 
described by St = (L/c)(t(11) — t(01)) An. As shown in Fig. 1, t is a 
function of g and V and the maximum value of t(11) — t(01) = T 
decreases as g is decreased from g = <x> to g = 2. However, for profiles 
with large index dips, T is constant with respect to g. This phenomenon 
is illustrated in Fig. 7e which plots T vs g for profiles with dips with 
various depths, d, and rms widths, w. Note that T ~ 0.155 for profiles 
with 100-percent deep dips (d = 1) whose widths are 25 percent of the 
core radius (w = 0.25). This explains why the maximum 5t in Fig. 6 is 
nearly the same for fibers S-l, G-l, and G-2, which had approximately 
the same A/i but different effective ^-profile exponents with deep and 
wide dips. 

Fibers G-l and G-2 were fabricated with intentional index grading, 
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Fig. 7 — The maximum delay lag of the lp(11) mode relative to the lp(01) mode 
plotted vs profile-^ for profiles with dips characterized by various depths, d, and 
widths, w. 
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having ^-profiles with successively smaller ^-values. Their delay dif- 
ference curves in Fig. 6 are horizontally displaced from the step-index 
curves and show that decreasing g increases the lp(01)-lp(11) equal- 
ization V-number and therefore the fiber diameter corresponding to 
zero intermodal dispersion. Fiber G-2 had an effective g ~ 3 and its 
equalization V = 4.25 number was significantly greater than the 
effective cut-off V co = 2.9 number for fabricated step-index fibers. 

The ( ) delay curves in Fig. 6 are computed by St = (L/ 

c)(t(11) - (t(01))Aai from the results in Fig. lb for fibers vrithg = 2.5 
power law profiles which are more optimal for double-mode operation. 
One of the curves applies to A = 0.011 core-cladding refractive-index 
differences corresponding to the described fabricated profiles. The 
other curve, with the smaller slope, applies to fibers with A = 0.001 
that are typical of currently fabricated large-core single-mode fibers. 
Intermodal dispersion tolerances as a function of fiber profile shape, 
V-number, source wavelength, and spectral bandwidth changes are 
related to the slope of the delay curves in the vicinity of their zero 
crossings. The curve in Fig. 6 for g = 2.5, V = 5.5, A = 0.001 indicates 
that 5t = 0.114 SV ns/Km. Changes in V-number are related to 
changes in wavelength through SV = -V [(1.28) /A 2 ] 8X. Therefore St 
= 490 8X ps/Km and a pulse broadening requirement of St = 1 ps/Km 
limits the source spectral bandwidth to S\ < 0.002 /zm which is typical 
of semiconductor lasers. From the curve in Fig. 3, for V = 5.5, g = 2.5, 
and no dip, Sg = —6.5 SV. Thus changes in profile g are related to 
changes in wavelength through Sg ~ 2.8 SX. Therefore a decrease of 
Sg « 0.1 would cause the zero intermodal dispersion wavelength to 
decrease by S\ ~ 0.036 /xm. This wavelength change could be compen- 
sated for by decreasing the fiber diameter and therefore the fiber V- 
number. 

IV. CONCLUSIONS 

We have reported the first time-domain transmission measurements 
in double-mode fibers which verify that intermodal dispersion effects 
can be very small in the 1.3-/un wavelength region where material 
dispersion effects are also very small. 

One of the graded-index fibers fabricated for this study, with an 
effective profile g ~ 3, had zero intermodal dispersion at A = 1.24 fim. 
This fiber had a V = 4.25, which is significantly greater than the 
effective cut-off V co = 2.9 for fabricated step-index single-mode fibers. 
Exact numerical solutions of the scalar wave equations were used to 
predict that fibers fabricated with g = 2.5 power law profiles without 
a central index dip have zero intermodal dispersion if the fiber 
F-number is 5.5. Additional numerical results showed that profile dips 
could be partially compensated for by reducing the profile g parameter. 
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For a system wavelength near X = 1.3 /xm, V = 5.5 corresponds to a 
0.1 -percent core-cladding index difference within a core diameter of 35 
/xm. This is approximately twice as large as the core diameter of a step- 
index single-mode fiber designed for the same wavelength." 
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